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Abstract 

We present the current status of the application of our approach of exact amplitude-based 
resummation in quantum field theory to precision QCD calculations, by realistic MC 
event generator methods, as needed for precision LHC physics. In this ongoing program 
of research, we discuss recent results as they relate to the interplay of the attendant IR- 
Improved DGLAP-CS theory of one of us and the precision of exact NLO matrix-element 
matched parton shower MC's in the Herwig6.5 environment in relation to recent LHC 
experimental observations. There continues to be reason for optimism in the attendant 
comparison of theory and experiment. 
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1 Introduction 



With the recent announcement [T] of an Englert-Brout-Higgs (EBH) |2J candidate boson 
after the start-up and successful running of the LHC for 2.5 years, the era of precision 
QCD, by which we mean predictions for QCD processes at the total precision tag of 1% or 
better, is squarely upon us. The attendant need for exact, amplitude-based resummation 
of large higher order effects is now more paramount, given the expected role of precision 
comparison between theory and experiment in determining the detailed properties of the 
newly discovered EBH boson candidate. Three of us(B.F.L.W,S.K.M,S.A.Y) have argued 
elsewhere [3111] that such resummation allows one to have better than 1% theoretical 
precision as a realistic goal in such comparisons, so that one can indeed distinguish new 
physics(NP) from higher order SM processes and can distinguish different models of new 
physics from one another as well. In what follows, we present the status of this approach to 
precision QCD for the LHC in connection with its attendant IR-improved DGLAP-CS [515] 
theory 0E] realization via HERWIRI1.031 [9] in the HERWIG6.5 [W\ environment in 
interplay with NLO exact, matrix element matched parton shower MC precision issues. 
We will employ the MC@NLO [UJ methodology to realize the attendant exact, NLO 
matrix element matched parton shower MC realizations for both HERWIRI1.031 and 
HERWIG6.5 in our corresponding comparisons with recent LHC data that we present 
herein. 

The discussion will therefore be seen to continue the strategy of building on exist- 
ing platforms to develop and realize a path toward precision QCD for the physics of the 
LHC. We exhibit explicitly a union of the new IR-improved DGLAP-CS theory and the 
MC@NLO realization of exact NLO matrix element (ME) matched parton shower MC the- 
ory. As our ultimate goal is a provable precision tag on our theoretical predictions, we note 
that we are also pursuing the implementation [T2] of the new IR-improved DGLAP-CS the- 
ory for HERWIG++ [13], HERWIRI++, for PYTHIA8 [Hj and for SHERPA PJ, as well 
as the corresponding NLO ME/parton shower matching realizations in the POWHEG [16] 
framework. For, one of the strongest cross checks on theoretical precision is the difference 
between two independent realizations of the attendant theoretical calculation. Such cross 
checks will appear elsewhere [T2] . 

In order to expose properly the interplay between the NLO ME matched parton shower 
MC precision and the new IR-improved DGLAP-CS theory, we set the stage in the next 
section by showing how the latter theory follows naturally in the effort to obtain a provable 
precision from our approach [4] to precision LHC physics. In the interest of completeness, 
we review this latter approach, which is an amplitude-based QEDcg)QCD(= QCDcg)QED) 
exact resummation theory [4] realized by MC methods, in the next section as well. We 
then turn in Section 3 to the applications to the recent data on single heavy gauge boson 
production at the LHC from the perspective of the analysis in Refs. [9] of the analogous 
processes at the Tevatron, where we will focus in this paper on the single Z/7* production 
and decay to lepton pairs for defmiteness. The other heavy gauge boson processes will be 
taken up elsewhere [T2]. Section 4 contains our summary remarks. 
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2 Recapitulation 



The starting point for what we discuss here may be taken as the fully differential repre- 
sentation 

da = / dxidx2Fi(xi)Fj(x2)da Tes (xiX2s) (1) 
i,3 J 

of a hard LHC scattering process using a standard notation so that the {Fj} and da res 
are the respective parton densities and reduced hard differential cross section where we 
indicate the that latter has been resummed for all large EW and QCD higher order 
corrections in a manner consistent with achieving a total precision tag of 1% or better 
for the total theoretical precision of ([1]). The key theoretical issue for precision QCD for 
the LHC is then the proof of the correctness of the value of the total theoretical precision 
Ac"th of (JTJ) . This precision can be represented as follows: 

Affth = AF © Aa res (2) 

where AA is the contribution of the uncertainty on the quantity A to Acr t J§. In or- 
der to validate the application of a given theoretical prediction to precision experimental 
observations, for the discussion of the signals and the backgrounds for both Standard 
Model(SM) and new physics (NP) studies, and more specifically for the overall normal- 
ization of the cross sections in such studies, the proof of the correctness of the value of 
the total theoretical precision Acx t h is essential. If a calculation with an unknown value of 
A<7 t h is used for the attendant studies, the NP can be missed. This point simply cannot 
be emphasized too much. 

In the interest of completeness here, we note that, by our definition, Acr t h is the to- 
tal theoretical uncertainty that comes from the physical precision contribution and the 
technical precision contribution [T7]: the physical precision contribution, Acr t p h hys , arises 
from such sources as missing graphs, approximations to graphs, truncations,....; the tech- 
nical precision contribution, A<7*{f h , arises from such sources as bugs in codes, numerical 
rounding errors, convergence issues, etc. The total theoretical error is then given by 

Aa th = A< ys © A< h . (3) 

The desired value for A<r t h, which depends on the specific requirements of the observa- 
tions, as a general rule, should fulfill Aa t h < fAa expt , where A<7 exp t is the respective 
experimental error and / < ~. This would assure that the theoretical uncertainty does 
not significantly adversely affect the analysis of the data for physics studies. 

In order to realize such precision in a provable way, we have developed the QCD (g) 
QED resummation theory in Refs. jl] for the reduced cross section in (Cp) and for the 
resummation of the evolution of the parton densities therein as well. In the interest of 
completeness and also because the theory in Refs. plj is not widely known, we recapitulate 
it here briefly. Specifically, for both the resummation of the reduced cross section and 
that of the evolution of the parton densities, the master formula may be identified as 



2 Here, we discuss the situation in which the two errors in ([2]) are independent for dcfmiteness; ([2]) has 
to be modified accordingly when they are not. 
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fa _ p SUM m (QCED) V-oo _J_ r rrn rf3fc Ji 
UO res — e Z^Ti,m=0 n!m! J Hj'i=l fc^ 

rr m d3fc/j 2 f _^_^y(pi+gi-P2-92-Efcj 1 -Efe' J2 )+-DQCED 

Pn,m(ki, • • • , k n ; k[, . . . , k' m ) -^jr , (4) 

where <i<7 res is either the reduced cross section da Tes or the differential rate associated to a 
DGLAP-CS p2,|6] kernel involved in the evolution of the {Fj} and where the new (YFS- 
style [IB]) non-Abelian residuals ^ n>m (ki, . . . , k n ; k[, . . . , k' m ) have n hard gluons and m 
hard photons and we show the final state with two hard final partons with momenta 
P2, <?2 specified for a generic 2f final state for definiteness. The infrared functions 
SUMir(QCED), -Dqced are defined in Refs. [U[7J|8]. This simultaneous resummation 
of QED and QCD large IR effects is exact. 

We stress that in Refs. [THS] the methods we employ for resummation of the QCD 
theory have been shown to be are fully consistent with the methods in Refs. [T91I20] . 
What is shown in Refs. [7HH] is that the methods in Refs. [HIEO] give approximations to 

our hard gluon residuals /3 n ; for, the methods in Refs. [T5P0]. unlike the master formula in 
fllj), are not exact results. Specifically, the threshold-resummation methods in Refs. [H 
using the result that, for any function f(z), 
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dzz n ' 1 f{z) 







<(-) max \f(z)\, 
n ze[o,i] 



drop non-singular contributions to the cross section at z — > 1 in resumming the logs in 
n-Mellin space. The SCET theory in Refs. [20] drops terms of O(X) at the level of the 
amplitude, where A = yJh./Q for a process with the hard scale Q with A ~ .3GeV so that, 
for Q ~ lOOGeV, we have A = 5.5%. From the known equivalence of the two approaches, 
the errors in the threshold resummation must be similar. Evidently, we can only use 
these approaches as a guide to our new non-Abelian residuals as we develop results for 
the sub-1% precision regime. 

We want to stress that, as it is explained in Refs. [3], the new non-Abelian residuals 
Pm jn allow rigorous shower/ME matching via their shower subtracted analogs: 

(3m,n ^ (3m,n (5) 

where the /3 m , n have had all effects in the showers associated to the {Fj} removed from 
them. This naturally brings us to the attendant evolution of the {Fj}; for, in order to 
have a strict control on the theoretical precision in (0Q), we need both the resummation of 
the reduced cross section and that of the latter evolution. 

When the QCD restriction of the formula in (j3J) is applied to the calculation of the 
kernels, Pab, m the DGLAP-CS theory itself, we get an improvement of the IR limit 
of these kernels, an IR-improved DGLAP-CS theory [TIE] in which large IR effects are 
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resummed for the kernels themselves. The resulting new resummed kernels, P^j? are given 
in Ref. [7H9] and are reproduced here for completeness: 



where the superscript "exp" indicates that the kernel has been resummed as predicted 
by Eq. (jlj) when it is restricted to QCD alone, where the YFS [18] infrared factor is 
given by Fyfs(o) = e~ CEa /T(l + a) where Ce is Euler's constant and where we refer the 
reader to Refs. [7J[8] for the detailed definitions of the respective resummation functions 
1a, Sa, /a, A = q, G □. Cf(Cg) is the quadratic Casimir invariant for the quark(gluon) 
color representation respectively. These new kernels yield a new resummed scheme for 
the parton density functions (PDF's) and the reduced cross section: 



with the same value for a in (pQ) with improved MC stability as discussed in Refs. [9] 
there is no need for an IR cut-off 'fco' parameter in the attendant parton shower MC based 
on the new kernels. We point-out that, while the degrees of freedom below the IR cut-offs 
in the usual showers are dropped in those showers, in the showers in HERWIRI1.031, as 
one can see from (jlj), these degrees of freedom are integrated over and included in the 
calculation in the process of generating the Gribov-Lipatov exponents 7,4 in (jSj). We note 
also that the new kernels agree with the usual kernels at 0(a s ) as the differences between 
them start in 0{a 2 s ). This means the the NLO matching formulas in the MC@NLO 
and POWHEG frameworks apply directly to the new kernels for exact NLO ME/shower 
matching. 

For completeness, we feature in Fig. |2] the basic physical idea underlying the new 
kernels as it was already discussed by Bloch and Nordsieck |23j: an accelerated charge 
generates a coherent state of very soft massless quanta of the respective gauge field so 
that one cannot know which of the infinity of possible states one has made in the splitting 
process q(l) —> q(l — z) + G <S> G\ ■ ■ • ® Gi, £ = 0, • • • ,00 illustrated in Fig. [2J The new 

kernels take this effect into account by resummimng the terms O ((ac s ln(-p) ln(l — z)) n ) 



3 The improvement in Eq. (j6j) should be distinguished from the resummation in parton density evolution 
for the "z — > 0" Regge regime - see for example Ref. 21,22 . This latter improvement must also be taken 
into account for precision LHC predictions. 
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Fj, a -»■ Fj, a' for 
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g -> g(l - 2) + G®Gi • • • ® G/, £ = 0, • • • , 
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Figure 1: Bloch-Nordsieck soft quanta for an accelerated charge. 



when z — > 1 is the IR limit. As one can see in ([7]) and (pQ), when the usual kernels are 
used these terms are generated order-by-order in the solution for the cross section o in 
([I]) and our resumming them enhances the convergence of the representation in ([1]) for a 
given order of exactness in the in input perturbative components therein. We now turn 
to the illustration of this last remark in the context of the comparison of NLO parton 
shower/matrix element matched predictions to recent LHC data. 



3 Interpaly of NLO Shower/ME Precision and IR- 
Improved DGLAP-CS Theory 



The new MC HERWIRI 1.031 [9 J gives the first realization of the new IR-improved kernels 
in the HERWIG6.5 [10] environment. Here, we compare it with HERWIG6.510, both 
with and without the MC@NLO pj] exact 0(a s ) correction to illustrate the interplay 
between the attendant precision in NLO ME matched parton shower MC's and the new 
IR-improvement for the kernels where we use the new LHC data for our baseline for the 
comparison. 

More precisely, for the single production at the LHC, we show in Fig. [2] in panel 
(a) the comparison between the MC predictions and the CMS rapidity data [25] and in 
panel (b) the analogous comparison with the ATLAS Pt data, where the rapidity data 
are the combined e + e~ — [l~ fi + results and the pt data are those for the bare e + e~ case, 
as these are the data that correspond to the theoretical framework of our simulations - we 
do not as yet have complete realization of all the corrections involved in the other ATLAS 
data in Ref. [26J. These results should be viewed from the perspective of our analysis in 
Ref. [9] of the FNAL data on the single Z/'-f* production in pp collisions at 1.96 TeV. 

Specifically, in Fig. 11 of the second paper in Ref. [9], we showed that, when the 
intrinsic rms parameter PTRMS is set to in HERWIG6.5, the simulations for 
MC@NLO/HERWIG6.510 give a good fit to the CDF rapidity distribution data [28] 
therein but they do not give a satisfactory fit to the DO distribution data [29] therein 
whereas the results for MC@NLO/HERWIRI1.031 give good fits to both sets of data 
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Figure 2: Comparison with LHC data: (a), CMS rapidity data on (Z/7*) produc- 
tion to e + e~, yU + /i~ pairs, the circular dots are the data, the green(blue) lines are 
HERWIG6.510(HERWIRI1.031); (b), ATLAS p T spectrum data on (Z/7*) produc- 
tion to (bare) e + e~ pairs, the circular dots are the data, the blue(green) lines are 
HERWIRI1.031(HERWIG6.510). In both (a) and (b) the blue(green) squares are 
MC@NLO/HERWIRI1.031(HERWIG6.510(PTRMS = 2.2GeV)). In (b), the green trian- 
gles are MC@NLO/HERWIG6.510(PTRMS =0). These are otherwise untuned theoretical 
results. 

with the PTRMS = 0. Here PTRMS corresponds to an intrinsic Gaussian distribution in 
Pt- The authors of HERWIG [27] have emphasized that to get good fits to both sets of 
data, one may set PTRMS = 2 GeV. Thus, in analyzing the new LHC data, we have set 
PTRMS = 2.2GeV in our HERWIG6.510 simulations while we continue to set PRTMS=0 
in our HERWIRI simulations. 

Turning now with this perspective to the results in Fig. [21 we see a confirmation of the 
finding of the HERWIG authors. To get a good fit to both the CMS rapidity data and the 
ATLAS p T data, one needs to set PTRMS = 2GeV [30] in the MC@NLO/HERWIG6510 
simulations. We again see that at LHC one gets a good fit to the data for both the rapidity 
and the p T spectra in the MC@NLO/HERWIRI1.031 simulations with PTRMS = 0. In 
quantitative terms, the x 2 /d.o.f. for the rapidity data and pr data are (.72, .72) ((.70,1.37)) 
for the MC@NLO/HERWIRI1.031(MC@NLO/HERWIG 6510(PTRMS=2.2GeV)) simu- 
lations. For the MC@NLO/HERWIG6510(PTRMS=0) simulations the corresponding 
results are (.70,2.23). 

Thus, we see that the usual DGLAP-CS kernels require the introduction of an intrinsic 
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Gaussian spread in p? inside the proton to reproduce the LHC data on the px distribution 
of the Z/'y* in the pp collisions whereas the IR-improved kernels give in fact a better fit 
to the data without the introduction of such a hard intrinsic component to the proton's 
constituents. This PTRMS is entirely ad hoc; it is in contradiction with the results of 
all successful models of the proton wave- function [3~I] . More importantly, it contradicts 
the known experimental observation of precocious Bjorken scaling [32j|33], where the 
SLAC-MIT experiments show that Bjorken scaling occurs already at Q 2 = 1 + GeV 2 for 
Q 2 = — q 2 with q the 4-momentum transfer from the electron to the proton in the famous 
deep inelastic electron-proton scattering process whereas, if the proton constituents really 
had a Gaussian intrinsic px distribution with PTRMS = 2GeV, these observations would 
not be possible. What can now say is that the ad hoc PTRMS = 2.2GeV value is really 
just a phenomenological representation of the more fundamental dynamics realized by the 
IR-improved DGLAP-CS theory. This raises the question of whether it is possible to tell 
the difference between the two representations of the data in Fig. |2j 

Physically, one expects that more detailed observations should be able to distinguish 
the two. Specifically, we show in Fig.EJthe MC@NLO/HERWIRI 1.031(blue squares) and 
MC@NLO/HERWIG6510(PTRMS=2.2GeV) (green squares) predictions for the Z/7* 
mass spectrum when the decay lepton pairs are required to satisfy the LHC type re- 
quirement that their transverse momenta {p T ,p T } exceed 20 GeV. We see that the high 
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Figure 3: Normalized vector boson mass spectrum at the LHC for p? (lepton) > 20 GeV. 

precision data such as the LHC ATLAS and CMS experiments will have (each already has 
over 5 x 10 6 lepton pairs) will allow one to distinguish between the two sets of theoretical 
predictions, as the peaks differ by 2.2% for example. Of course, other such detailed ob- 
servations may also reveal the differences between the two descriptions of parton shower 
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physics and we will pursue these elsewhere [12]. We await the release of the entire data 
sets from ATLAS and CMS. 

4 Conclusions 

We have shown that the realization of IR-improved DGLAP-CS theory in HERWIRI1.031, 
when used in the MC@NLO/HERWIRI1.031 exact 0{a s ) ME matched parton shower 
framework, affords one the opportunity to explain, on an event-by-event basis, both the 
rapidity and the pt spectra of the Z/^* in pp collisions in the recent LHC data from CMS 
and ATLAS, respectively, without the need of an ad hoc intrinsic Gaussian distribu- 
tion with rms value of PTRMS = 2 GeV in the proton's wave function. We argue that 
this can be interpreted as providing a rigorous basis for the phenomenological correctness 
of such ad hoc distributions insofar as describing these data using the usual unimproved 
DGLAP-CS showers is concerned and we have proposed that other distributions such as 
the invariant mass distribution with the appropriate cuts be used to differentiate between 
the fundamental description of the parton shower physics in MC@NLO/HERWIRI1.031 
and these phenomenological representations in MC@NLO/HERWIG6510. We have em- 
phasized that the precociousness of Bjorken scaling argues against the fundamental cor- 
rectness of the ad hoc ansatz with PTRMS = 2 GeV, as do the successful models [31] 
of the proton's wave function. We have the added bonus that the fundamental descrip- 
tion in MC@NLO/HERWIRI1.031 can be systematically improved to the NNLO parton 
shower/ME matched level which we anticipate is a key ingredient in achieving the sub- 
1% precision tag for such processes as single heavy gauge boson production at the LHC. 
Evidently, the use of ad hoc models would compromise any discussion of the theoreti- 
cal precision relative to what one could achieve from the fundamental representation of 
the corresponding physics via IR-improved DGLAP-CS theory as it is realized in HER- 
WIRI1.031 when employed in MC@NLO/HERWIRI1.031 simulations. We are pursuing 
additional cross checks of the latter simulations against the LHC data. 

In closing, two of us (A.M. and B.F.L.W.) thank Prof. Ignatios Antoniadis for the 
support and kind hospitality of the CERN TH Unit while part of this work was completed. 
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